Zero-energy vortex bound states in noncentrosymmetric superconductors 
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We consider bound states at the vortex core of a non-centrosymmetric superconductor. We siiow 
that, despite the mixing of singlet and triplet order parameters, zero energy states survive within 
certain parameter space as in vortices of some chiral p-wave states. 
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Noncentrosymmetric superconductors, such as 
CePtaSi, Li2Pt3B, and MgioIrigBig, have drawn a great 
deal of interest in the past few years [IHi. Microscopi- 
cally, the degeneracy among the two pseudospin states 
at a given momentum, connected by the composite oper- 
ations of time-reversal and inversion transformations, is 
lifted since the latter symmetry is broken. Hence, Cooper 
pairs can no longer be classified in terms of singlets 
versus triplets, nor s-wave versus p-wave In CePtaSi, 
for example, the superconducting state is believed to 
have a triplet order parameter d(k) oc (fc^y — fcj/x) in 
addition to a singlet one It has been proposed that 
this mixing is responsible for some of the seemingly 
contradictory behaviors of these superconductors 0, Q. 
Moreover, some intriguing broken symmetry properties 
of these superconductors have been predicted iQj . 

Recently, bound state 0] within the odd winding num- 
ber vortex of the p-wave k^+iky superfluid has been a hot 
topic since among them there is a localized zero-energy 
state[l- 11|, which, in terms of its associated creation op- 
erator, is a self-hermitian Majorana fermion [l2j . Propos- 
als that utilize such zero-energy states for implementing 
the topological quantum computationjjl. [lo| are based on 
its unique properties including robustness against per- 
turbations from deformations of order parameters and 
nonmagnetic impurities Ill4l5l| . 

It is well known that the bound state spectrum asso- 
ciated with an isolated vortex with winding number one 
in an s-wave superconductor is E = [n + with the 
energy scale w of order A^/i?i?[7| (here n is an integer, 
A the energy gap and Ep the Fermi energy), while the 
spectrum in the A phase of p-wave superfluid '^He is given 
hy E — nu which includes a state at zero energy [8|. In 
this paper, we shall study the vortex bound states in non- 
centrosymmetric superconductors described by a mixture 
of singlet and triplet order parameters, which is inter- 
mediate between the previous two examples. First, we 
demonstrate that there are two such zero-energy states 
(for each p^, to be defined later) corresponding to the 
pure p-wave order parameter with d(k) oc (fc^y — fcj,x) 
presenting the combined up and down equal-spin pair- 
ings. Thus the two states differ by spin orientation, and 
at first sight, would be coupled by the additional s-wave 
order parameter and hence acquire finite energies. How- 
ever, we shall show that, as long as the s-wave component 



is smaller than a critical value, the zero-energy states sur- 
vive. Besides, we note that this condition is identical to 
the existence of a nodal gap. In addition, we consider a 
Rashba spin-orbital interaction. We shall show that the 
zero-energy states again survive. 

The creation operator for a quasiparticle excitations 
in an inhomogeneous superconductor is a linear combina- 
tion of electronic annihilation and creation field operators 
and tpl^^{r)- 

("t«."im)(;j;y+(..m.n(n)(*;«) 

(1) 

satisfying [i/cft,a^] — ea^ where i?cff is the effective 
mean-field Hamiltonian and e is the quasiparticle energy. 
The coefficients {u,v) = (u-j-, m^, w^) satisfy the Bo- 
goliubov deGennes (BdG) equation, 



H 



(2) 



where, for the ordinary cases, Hq = — ^ — Ep is the 
kinetic energy (we shall add the possible Rashba interac- 
tion later). H is a two- by- two matrix due to the pairing. 
In the singlet case H represents As{r){ia2)- In the triplet 
case H = i(V-z3) + i5- V. The vector D = -iVkA(f,k) 
where A(r, k) is the order parameter parameter. By the 
conventional notation, A(r, k) = Ap(r)d(k) ■<j{ia2)- H is 
then a sum of the above two when both order parameters 
are present. 

We first consider the simpler case (c./. e.g. [l3|) where 
the order parameter is of pure p-wave character with 
d(k) = {ko,y - kyx)/pF. (Here pF = {2mEpy/^). We 
shall show that, for each value of momentum along the 
vortex line Pz less than pp, there are two zero-energy 
states with the associated wavefunction (u(r), {)(r))^ 
given by 



(l,0,-l,0)^i^l(p)e^P^^ 
(0,e*^0,-e-''^)^i?2(p)e*P^" , 



(3) 



in cylindrical coordinates f— {p,(j),z), where the radial 
functions i?i^2 are independent, finite, and decaying at 
infinity. 
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For an isolated vortex line with winding number 1 , the 
order parameter can be expressed as Ap(r) = Ap(p)e"^ 
where we shall choose the gauge where Ap(p) is real and 
positive. The coupling 11 in © is then 



Ap(p) 

PF 



(4) 

Ap(p) is zero at p = and increases toward its 
asymptotic value Aq within a range of coherence length 
? ^ Pp^- In principle we should add also respectively 
■^dpA{p) and ■^e^^'^dpA{p) to the upper left and 

lower right elements in eq Q due to the V • D term 
in n, but it can be shown that, since these terms are reg- 
ular p and vanish necessarily a.s p oo, they do not 
affect the arguments below and so we would not show 
them explicitly to simplify the equations. It is clear from 
^ that the BdG equation ([2]) for excitations with up 
and down spins become decoupled. For the states with 
e = 0, we denote the two independent excitations by 
and aj. For a\ with up spin, the wavefunctions can be 
factored into u^ir) = e'^P'^Uf{p) and v-f{r) = e^^'^v^{p). 
The equations for Uf + and — Vf become decoupled: 



2m ^ dp^ 



Id, A 
p dp PF 



pjp) , d 
^dp 



(5) 
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where fj.{pz) = Ep — |^ is an effective chemical poten- 
tial. Here we assume Pz < Pfi a-nd thus p > 0. Two 
properties about Eq. ([S|) have to be noted here. First, 
for each of the above equation the presence of a regular 
singular point at the origin forces the general solutions 
to be divergent at the origin. Second, the equation cor- 
responding to u-[ + v-[ with a minus sign in front of d/dp 
has its both solutions unbounded at infinity, whereas the 
equation for u^ — Vf has two solutions decaying as p — oo. 
The above can be shown by substituting e'^n'' into the 
asymptotic form of (O where terms proportional to l/p 
can be neglected and Ap(p) is replaced by Ap. Then 
we have —x"^ + -gj =F = 0, where x = p\\/pf- For 

the lower sign appropriate to u-f — Wf , the solutions for x 

and hence p|| is given by PF[+i^^ ± (assuming a 

weak-coupling superconductor and thus Aq <C Ep) with 
positive imaginary parts (if /2 > 0) . Hence u-f — w-f- has a 
pair of decaying solutions xi(p) ^^^d x^ip)- For the up- 
per sign appropriate to u-\- + v-f, we only have two growing 
solutions that must be rejected and so u-f{p) + v^[p) = 
must hold. Now both solutions Xi(p) and x^ip) in gen- 
eral contain a divergence of In p as p —J' 0. Nevertheless, 
a finite solution as p — can be found by an appropriate 
linear combination so that we can write the zero-energy 
state as the first line in ([3]), -Ri(p) = axi{p) + ^X2(p) 
where the coefficients are determined to cancel In p near 



the origin. (For \pz\ > pf, A < 0, then both ± 
have a single solution which decays at infinity and thus 
no zero-energy state is allowed, c./., e.g., lol 1^). 



The state associated with a\ can be obtained in a 
similar manner. The factorization is {ui{r) , ui{r)) = 
e'^P'^ {e^'^ui{p) , e~^'^vi{p)). The equations are also decou- 
pled for uizLvi as in previous case, but the kinetic energy 
contains an additional term — 1/p^, which causes a diver- 
gence of 1/p for the solution near the origin. Analogous 
procedure shows that there are two decaying solutions 
771 (p) and 772 (p) for — -D^, but two exponentially in- 
creasing solutions for ui + v^. Therefore we have the 
solution as (u2(r), ^2(r)) = 6*^-^(0, e*"^, 0, -e-*'^)i?2(p), 
where R2 = crji + dr]2 is a suitable linear combination to 
cancel the divergence of 1/p. Now we obtain the eigen- 
functions associated with the two zero-energy excitations. 
We note the relation 



U2 — e"^(TiUiK(p) , V2 = e "^(TiWik(p), 



(6) 



where k(p) = R2{p)/Ri{p), which will be useful later. 

With the solutions ^ obtained for a pure p-wave or- 
der parameter ([4]), we are going to consider the effects of 
lacking inversion symmetry on ([3]) . For simplicity of pre- 
sentation, we shall consider separately the Rashba inter- 
action and an admixture of singlet order parameter but 
we only state the general conclusion at the end. First, we 
shall show via perturbation theory that a small Rashba 
interaction or a small s-wave order parameter would not 
destroy the zero energy states. Then we shall consider 
general magnitude of these two interactions. 

Now we add a Rashba spin-orbital interaction h = 
{—az X p • (?) to the the kinetic energy parts Hq in ([2]). 
In cylindrical coordinates, h and its counterpart —h* as- 
sociated with the hole sector can be written as 







(7) 



where the upper (lower) sign is for the electron(hole) sec- 
tors, respectively. The expectation values of the this 
spin-orbital interaction for either of the two states given 
above are obviously zero. The matrix element between 
the two states in Q is proportional to the spatial in- 
tegral of u\hui — v\h*vi, which, with eq ([B]), equals 
{u\ai[e^^'f'h -f e"^(— ft.*)]wi} times some function of p. 
Hence the matrix element is zero by explicit use of ([7|. 
Therefore, within perturbation, the two states in ^ are 
unaffected by the coupling from the spin-orbital interac- 
tion. 

Now consider an additional a singlet pairing order pa- 
rameter. We therefore add a term 



A.(r') = 



e"t' 
-e^-^ 



A.(p) 



(8) 
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to the pure triplet one in (|4]) in eq ([2]) . For mixing which 
do not break time reversal symmetry far from the vor- 
tex, the ratio limp_j.oo {As{p) / Ap{p)) = /3 must be real. 
The expectation value of ([5]) in any given state in ^ 
is again obviously zero, and matrix element of (|5|) be- 
tween the two states in ^ involves the spatial integral 
oi u\e^'^ {ia2)vi -f U2e~*'^(— 10-2)"! which, by using ([6|), is 
zero. Hence in the small /3 regime, the two zero-energy 
states remain. 

The previous two paragraphs demonstrate that the two 
states in ([3]) remain according to perturbation theory. 
Now we consider the more general case. We shall show 
that there are two zero-energy states with general form 



(9) 



with u-j- = —v^ and uj, = —v^ all finite and decaying at 
infinity, which survive the additional interactions. 

First consider order parameter mixing. The operator 
n in BdG equation is a sum of (jlj and The corre- 
sponding set of differential equations are: 



Lqu-^ - Put ~ ^sip)vi = , 
Lqv^ - Pu^ - A.s{p)ui = , 

iin-P% + A.(p)ut = 0, (10) 



where the differential operators L„ — -^{-^ 



■ p, and P 



A(p) , 
PF ^ dp 



\_d_ 

P dp 



(dB + 2^)- Denoting = ± 



and — ± v^, the above equations for the 's are 
decoupled with the w~ 's. Writing 
= iwj^ , we arrive at 



w^± iw^ and 



(Lo + P- 
(Lo-P 



1 



T iAs{p))z^ - 



1 



1 



Amp 
1 



Amp- 



(11) 

,2^^T=0- (12) 



Though in general (z^) couples with zZ (2I), we 
note that at infinity the above set of equations for each 
of and z^ couples to itself only. Then for a given 
Pz — Pf cos 9 associated with the excitation, we can write 
p = Ep sin^ > 0. satisfies, asymptotically for large 



are given by PF[i-:^^ i y sin^ ± The imaginary 

parts are therefore, for small Ao/Ep, ^-§^(1 ± ^j^)- 
For a given and 6, there are four roots associated with 
z^ with positive imaginary parts when 4^ < 1, which 
leads to four independent decaying solutions for z^ . The 
same arguments show that the solutions for z^ are expo- 
nentially increasing as p — >■ cx), and hence we must choose 
— — 0, that is, Uf = —Vf and uj, = —5^ as in ([9|). 
Therefore a general solution for the zero-energy state for 
Eq. (jlO[) can be represented by the linear combination of 
the 4 decaying solutions. 



/ Mp) \ 
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uiip) 
v^{p) 




yi{p) 
-1 


V hip) J 


1=1 


^ -y.[p) 1 



Mp) 



(14) 



where the /'s decays towards zero at infinity, and so are 
the Uifi. From eq PH]) . the divergences near the origin 
are of the form In p for the first row, and 1 /p for the sec- 
ond row. Now we have two equations determining the c's 
by which the respective divergence can be removed. In 
general, we can have two independent sets of {ci} satis- 
fying the above, which in turn leads to two independent 
zero-energy states within the vortex core. A crucial con- 
sequence is drawn from the above arguments. For the 
relative pairing strength |/3| > 1, |/3/sin6'| > 1 so that 
the zero-energy states no longer survive, which we con- 
clude that, in addition to /i = 0, |/?c| = 1 is another 
critical parameter. On the other hand, the zero-energy 
states exist at the core when |/3| < 1. The density of 
such excitations (number per unit length of the vortex 
line) can be determined from the condition |/3/ sin0| < 1. 
For a spherical Fermi surface, we obtain ^f^y^l — /?^. 
We note that the energy gap for the bulk excitations are 
given by Ao|sin6' ± /3|, hence the above critical value of 
6 for the existence of the E = vortex bound state cor- 
responds to exactly the existence of a nodal line in one 
of the branches. This is reasonable as this is the value 
of Pz where two of the decaying solutions for s'mO > |/3| 
become extended, destroying the possibility of obtaining 
the solution eq © which converges both at p — > 00 and 
0. 

Next we move on to the case when the spin-orbital in- 
teraction ([7]) is included in ^ . With again the wavefunc- 
tion in the form of eq ([9]) , the zero-energy BdG equation 
can be written as 



^^ + ^^+EpsinHT^PA, 
2m dp Pf dp 



(13) 



where one recalls that As{p — ?> 00) = /3Ao. Note that a 
similar equation as ()13p except the positive coefficient in 
front of d/dp is for the z+'s. Again, take e*Pii'' as asymp- 
totic solutions. In the weak-coupling limit, pu for z~ 



Lou-f - Put ^ "(-7- + -)h = : 
dp p 

Lqv^ ~ Pu^ - a{— + -)vi = , 
dp p 

d 

Liui - Pvi + a—u^ = , 
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Pu, + a-v, 







(15) 



We can analyze these equations in the same manner as 
the previous case. At infinity the BdG equations become 
decoupled as 



1 d 



2m dp 



PF 



dp 



(16) 



where cos 9 — p^/pF, and note there is a similar equation 
for except the overall positive coefficient associated 
with d/dp. p\\'s associated with the asymptotic solution 
e'PiiP have imaginary part ^^^(Izfc . "^'"^ ), which 

is positive for all 9. (Here vp = pp/m). On the contrary, 
the corresponding p|| for have only negative imaginary 
parts. Applying the same arguments as before, the zero- 
energy bound states survive under any magnitude of the 
spin-orbital interaction. One can understand this result 
by the fact that the size of the Fermi surfaces at Pz for the 

1/2 

two branches are given hy pp± = ^{2mp) + (m^a^)] ± 

ma = pf \J i^)"^ + sin^ 9 ± which remain finite for 
arbitrary large value of a. 

The above analysis can be generalized to the case 
where both a and /3 are finite. We find that a pair of 

E 



states exist if one has both ^ + ^in^ 9 + 

! > 0. We note 



-h /3 > and ./( — )2 -t-sin2 6' - ^ 

that in the helicity basis (spin quantization axis along 
{z X p)), the order parameter on the ± branches of the 
Fermi surfaces are given respectively by Ap^^^ ± As = 



hence the existence or 

absence of the E = Q bound states is determined by the 
relative sign of the order parameter on these two Fermi 
surfaces. 

The solutions in ([9]) taking the lack of inversion into 
account are as robust as that of ^ in p-wave superfluids. 
It is evident that the local charge density is zero since the 
solutions have equal magnitudes for electron and hole ex- 
citations of the same spin projection. Hence the states 
are not susceptible to nonmagnetic impurities. Similarly, 
the states are not affected by the Zeeman magnetic field 
along z again due to the particle-hole symmetry. Further- 
more, they are also not altered by exchange or Zeeman 
fields in the in-plane directions, for the azimuthal depen- 
dence in each of ([9]) leads to zero matrix elements of the 
local spin density operator among the states [l5]. 

With the two general independent solutions of the 
form in the corresponding creation operators are not 
necessarily self-hermitian. Here we show that a set of 
two independent Majorana fermions can be built from 
them. We first demonstrate this for p^ =0. From 
the two linearly independent solutions, one first con- 
struct two orthonormal wavevectors (u, v) and thus two 



corresponding operators a[ and a\. From orthonormal 
properties one then have {a|,ai} = {a2,a2} = 1 and 
{a|,a2} — 0. Since [Hoff,ai] = Eia\ = 0, we also have 
[-ffcffjCKj] — —EiOi — for i = 1,2. Since we have only 
two linearly independent solutions to eq ([2]) of zero ener- 
gies, a\ 2 must just be a linear combinations of cki,2. We 
denote this in matrix notation as ~ Ca where C is a 
2x2 matrix. We shall find the transformation 7^ = Wa^ 
such that the 7's are independent Majorana fermion op- 
erators, that is, 7,| — 7i for i — 1,2, and {71, 7]} = 0. By 
choosing W as unitary, we can assure the conditions of 
normalization and orthogonality for the 7's. It remains 
only to make the 7's self-conjugate. From {ai, a]} = ^ij 
where i,j = 1 or 2, one can show that C is unitary. 
Furthermore, — CC*a^ and thus C^^ = C*. Since 
C is unitary, C is also symmetric. Hence we can write 
C = e'''*e'"" '^, where A is real and h is perpendicular to 
y. We thus have 7^^ = WCa = WCW"^7 which would 
equal 7 if we choose W = g^'f "■'^e^*'^/^^ Thus the op- 
erators 7^ , 72 constitutes a set of two independent Majo- 
rana fermions. The constructions of Majorana fermions 
for finite pz^s can proceed in a similar manner if we re- 
place the e'P^^ factors in the wavef unctions by cos{pzz) 
and si'a{pzz). 

As demonstrated already in, e.g., the dimension 

of Hilbert space of a Majorana fermion is \/2, that is, each 
two Majorana fermions combines to form one Fermionic 
state with two degrees of freedom (occupied or empty). 
Hence, for our system with n„ vortices per unit area, we 

have a residual entropy density \Jl — P"^ ~ In 2 . 
This ground state degeneracy is lifted only by the finite 
overlap between the vortices [l4| . the resulting energies 
are thus exponentially small in the vortex spacings. The 
existence of this residual entropy can be used to demon- 
strate the existence oiE = Q vortex bound states, as well 
as be a measure of the mixing of the two superconducting 
order parameters. 

In conclusion, we have considered the vortex bound 
states in a non-centrosymmetric superconductors, in par- 
ticular for an order parameter appropriate to CePtaSi. 
We demonstrated that the zero energy states exist only 
for certain range of pz values depending on the magni- 
tudes of the singlet versus the triplet order parameters. 

While preparing this manuscript, we become aware of 
17| which however discusses a very different aspect of 
vortex bound states of non-centrosymmetric supercon- 
ductors. 
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